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The Bogoliubov approximation is used to study the excited states of a 
dilute gas of N atomic bosons trapped in an isotropic harmonic potential 
characterized by a frequency ujq and an oscillator length do = \fKJmuo~. The 
self-consistent static Bose condensate has macroscopic occupation number 
JVo > 1, with nonuniform spherical condensate density no(r); by assump- 
tion, the depletion of the condensate is small {N' = N — Nq <C Nq). The 
linearized density fluctuation operator ft and velocity potential operator <£' 
satisfy coupled equations that embody particle conservation and Bernoulli's 
theorem. For each angular momentum I, introduction of quasiparticle oper- 
ators yields coupled eigenvalue equations for the excited states; they can be 
expressed either in terms of Bogoliubov coherence amplitudes ui(r) and vi{r) 
that determine the appropriate linear combinations of particle operators, or 
in terms of hydrodynamic amplitudes p'i(r) and 3>5(r). The hydrodynamic 
picture suggests a simple variational approximation for I > that provides 
an upper bound for the lowest eigenvalue lo\ and an estimate for the corre- 
sponding zero-temperature occupation number N[\ both expressions closely 
resemble those for a uniform bulk Bose condensate. 
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I. INTRODUCTION 



Recent experimental verification of Bose condensation in dilute confined 87 Rb |IJ has 
stimulated extensive theoretical research. In the Bogoliubov approximation, N atoms 
occupy the macroscopic condensate, with Nq < N and only small depletion (N'/N = 
1 — N Q /N <C 1). For a repulsive scattering length a > and mean atomic density n, 
this depletion is of order y/na 3 B and thus nonperturbative. 

The balance between kinetic energy h 2 /2m^ 2 and interaction energy Anah 2 n/m defines 
the "healing length" £ = (87ran) _1//2 , and the Bogoliubov approximation requires that a <C £. 
In addition, a confined Bose gas with typical dimension R differs qualitatively from an ideal 
Bose gas whenever £ <C R. The present work uses the Bogoliubov approximation to provide 
a hydrodynamic description of the excited states and a variational estimate of the eigenvalue 
u>i and zero-temperature occupation number N[ for the lowest radial mode for each positive 
angular momentum I > 1. 

II. BASIC FORMALISM 

In the harmonic trap U = \muj\r 1 (here taken as isotropic) has a size do = \Jh/ mu « 
1 /iin, and the positive scattering length a m 10 nm acts to expand the atomic cloud. For 
M = Na/do ^> 1, the actual atomic density ~ N/R 3 is smaller than the naive estimate 
~ N/d 3 ), with the radial expansion factor R/d of order A/" 1 / 5 0; thus the Bogoliubov 
approximation holds for Na 3 /R 3 w Af 2/5 (a/d ) 2 < 1 [namely, iV < (d /a) 6 ~ 10 12 ]. 
The additional condition R ^> ^ for the failure of an ideal-Bose-gas model now requires 
(SttA/" 4 / 5 ) 1 / 2 > 1, which here holds even for N w 100. 

The condensate wave function \l/ satisfies the nonlinear Gross-Pitaevskii equation [||,[5| 

(T+C/-/i + y)* = 0, (1) 

where T = — fr 2 V 2 /2m is the kinetic energy, V = 47ra/i 2 m~ 1 |\l/| 2 characterizes the repulsive 
interaction energy, and the chemical potential [i must be adjusted to ensure that / dV |^| 2 = 
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Nq and that ^ — > at large distances. In the present case, the condensate is stationary 
with density n (r) = |\I/(r)| 2 . For large M , the kinetic energy can be ignored ||, giving the 
"Thomas-Fermi" approximation 

V TF = {H~U) 0{h -U) = \huo {R 2 - r 2 ) 9{R - r), (2) 

where lengths are expressed in units of do, p = ^TlujqR 2 defines the dimensionless "radius" 
R of the condensate, and R 5 = 15JV ||. 



III. HYDRODYNAMIC DESCRIPTION 

The (small) noncondensate field operator </> obeys the linear equation 

ihd^/dt = (T + U - p + 2V) 4> + V<f>\ (3) 

along with the adjoint equation. In the Bogoliubov approximation (N' <C Nq), the total 
particle density operator becomes ffiip ~ uq + p' , where the operator p' = ^Jno{4> + $) 
characterizes the fluctuating noncondensate density; similarly, the particle current operator 
here becomes j' ~ n V$', where $' = (h/2miy/n^)(<j) — ft) is the velocity potential operator. 
Equation (||) readily yields 

dp'/dt + V • (n V$') = 0, (4) 



and the linearized Bernoulli's equation | 



+ 5-ff + 



v-(p'^)-W 

n 



0. (5) 



dt m 2 4m% 

The latter can be recognized as a linearization of the exact classical Bernoulli's theorem |7j 
U + \ttlv 2 + (e + p) /n + md<& /dt = for a compressible irrotational isentropic fluid with 
density n + p' and velocity potential $' [here, p(n) = 2iiah 2 n 2 jm is the pressure functional 
and e(n) = \JnT\fn + p(n) is the (quantum) energy-density functional]. 
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IV. BOGOLIUBOV EQUATIONS 



The hydro dynamic description provides a valuable qualitative picture of the normal 
modes of the nonuniform compressible spherical condensate, but the resulting equations 
involve spatial derivatives of no; furthermore, in the Thomas- Fermi approximation the 
appropriate boundary conditions at R are not obvious. Thus it is preferable to return to 
Eq. (d), expressing the field operator as a linear combination of quasiparticle operators <x,- 
and a] 

0(r, t) = [ujtfaje-** - «;(r)a}e***] ; (6) 
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here, the wave-function coefficients Uj and Vj constitute a two-component vector function 
Uj that obeys the matrix "Bogoliubov" eigenvalue equation || 

(T + U - u + 2V) Uj - V n Uj = hujj r 3 Uj. (7) 

In contrast to Eqs. @) and (||), only V oc n appears here, and the quantum interpreta- 
tion provides the obvious boundary condition that Uj vanish for r — > oo (note that the 
hydrodynamic amplitudes are essentially linear combinations of Uj and Vj). 

In the case of a static spherical condensate, the Bogoliubov Eqs. (|7|) have solutions of the 
form ui(r)Yi m and vi(r)Yi m ; in addition, they have a simple variational basis |3| that yields 
an upper bound for the lowest eigenvalue uji for each positive angular momentum I > 0. To 
ensure the proper normalization, take W/(r) = cosh xi fi{r) and Vi(r) = sinh xi fi{r), where 
Jq°° r 2 dr \fi{r)\ 2 = 1. With the Thomas- Fermi approximation (fj), the analogy to acoustic 
waves in a sphere suggests the radial trial function fi(r) oc x (1 — x 2 )8(l — x), where 



x = r/R. A simple variational calculation yields the estimate oui/ujq = yT 2 + 27] VJ, with 
Ti = | (2/ + 3) (21 + 7)/R 2 and Vi = 3R 2 / (21 + 9) the expectation values of the dimensionless 
kinetic energy — |V 2 and Thomas-Fermi potential energy \(R 2 — r 2 ). This frequency has a 
very different form for small and large / 



Ul | y 5 (2/ 4- 3) (2/ 4- 7) /(2/ 4- 9), for T t < Vf, 

( 8 ) 

Uo { (2Z + 3)(2/ + 7)/4R 2 , forT z >^. 
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Stringari || has obtained a related expression uji/ujq = \/l with a purely hydrodynamic 
model, adding that W\ = uj is an exact result. Note the close analogy to the Bogoliubov 
excitation energy = \Jt% + 2TkVk for a uniform bulk medium, where T k = h 2 k 2 /2m and 
V k = 47rah 2 n /m @. 

Each of these lowest excited states for I > has a zero-temperature occupa- 
tion number N[ = r 2 dr |t>;(r)| 2 , with the variational estimate N[ = sinh 2 xi — 
|(Tj + Vz)/ \Jt 2 + 2T;VJ — |, again very similar in structure to that for a uniform Bose 
gas §. In particular, iV/ w y/Vi/8Ti > 1 for T z < VJ, and AT/ w VJ 2 /4T 2 < 1 for 7} > VJ. 
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